Searching for non-Gaussianity in the WMAP data 
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Some analyses of recent cosmic microwave background (CMB) data have provided hints that 
there are deviations from Gaussianity in the WMAP CMB temperature fluctuations. Given the far 
reaching consequences of such a non-Gaussianity for our understanding of the physics of the early 
universe, it is important to employ alternative indicators in order to determine whether the reported 
non-Gaussianity is of cosmological origin, and/or extract further information that may be helpful 
for identifying its causes. We propose two new non-Gaussianity indicators, based on skewness and 
kurtosis of large-angle patches of CMB maps, which provide a measure of departure from Gaussianity 
on large angular scales. A distinctive feature of these indicators is that they provide sky maps of 
non-Gaussianity of the CMB temperature data, thus allowing a possible additional window into 
their origins. Using these indicators, we find no significant deviation from Gaussianity in the three 
and five-year WMAP ILC map with KQ75 mask, while the ILC unmasked map exhibit deviation 
from Gaussianity, quantifying therefore the WMAP team recommendation to employ the new mask 
KQ75 for tests of Gaussianity. We also use our indicators to test for Gaussianity the single frequency 
foreground unremoved WMAP three and five-year maps, and show that the K and Ka maps exhibit 
clear indication of deviation from Gaussianity even with the KQ75 mask. We show that our findings 
are robust with respect to the details of the method. 

PACS numbers: 98.80.Es, 98.70.Vc, 98.80.-k 



I. INTRODUCTION 

Within the standard approach to cosmological mod- 
elling the suggestion that the Universe underwent a 
brief period of rapid acceleration expansion [1] before 
the epoch of primordial nucleosynthesis has become an 
essential building block of the standard cosmological 
model. Besides solving the so-called flatness, horizon and 
monopole problems, such a period of cosmological infla- 
tion provides a mechanism for the production of the pri- 
mordial density fluctuations, which seeded the observed 
cosmic microwave background (CMB) anisotropics and 
the formation of large-scale structure in the Universe. 

There are more than one hundred inflationary mod- 
els (see, e.g., the review articles Refs. [2]), among which 
the simple ones are based on a slowly-rolling single 
scalar field. An important prediction of a number of 
these simple models is that they can generate only tiny 
non-Gaussianity, which should be undetectable in the 
Wilkinson Microwave Anisotropy Probe (WMAP) CMB 
data 0. There are, however, a large class of inflationary 
models that can generate non-Gaussianity at a level de- 
tectable by WMAP [4] . These scenarios comprise models 
based upon a wide range of mechanisms, including special 
features of the inflation potential, multiple scalar fields, 
non-canonical kinetic terms, and non-adiabatic fluctu- 
ations (see the review Ref. ^ and references therein). 
Thus, the detection of non-Gaussianity in CMB data may 
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potentially be useful to discriminate inflationary models 
and shed light on the physics of the early universe. 

In the statistical analyses by using one, three and five- 
year 0, 0, 0] CMB measurements along with some dif- 
ferent statistical tools, the WMAP team have found that 
the CMB data are consistent with Gaussianity. However, 
some recent analyses have provided clear hints that there 
are significant deviations from Gaussianity in the WMAP 
data. Clearly the study of detectable non-Gaussianities 
in the WMAP data must take into account that they may 
have non-cosmological origins as, for example, unsub- 
tracted contamination from galactic diffuse foreground 
emission [lo', Tl] and unconsidered point sources |121]. If 
they turn out to have a cosmological origin, however, this 
could have far-reaching consequences on our description 
of the Universe, particularly on the inflationary picture. 

In view of this, a great deal of effort has recently gone 
into verifying the existence of such non-Gaussianity by 
employing several different statistical signatures of non- 
Gaussianity in its various forms (see, e.g., Refs. [l^ and 
related Refs. [l3])- Apart from revealing the existence of 
non-Gaussianity in CMB data, different statistical tools 
are sensitive to different systematics and may be useful 
in determining their origins. In addition, different in- 
dicators can in principle provide information about the 
multiple types of non-Gaussianity that may be present in 
CMB data. It is therefore important to test the data for 
deviations from Gaussianity by using a range of different 
statistical tools to identify any non-Gaussian signals on 
the CMB sky. 

In this paper, we propose new large- angle non- 
Gaussianity indicators, based on skewness and kurtosis 
of large-angle spherical-shaped patches of CMB maps. 
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which provide a measure of departure from Gaussianity 
on large angular scales. A distinctive feature of these 
indicators is that they provide sky maps (directional in- 
formation) of non-Gaussianity of the CMB temperature 
fluctuations, thus allowing a possible additional window 
into their causes. Using these indicators, we find no sig- 
nificant deviation from Gaussianity in the WMAP three 
and five-year ILC KQ75 masked map, but ILC unmasked 
map exhibit deviation from Gaussianity. On the other 
hand, our indicators reveal deviations from Gaussianity 
of different degree in the three and five-year single fre- 
quency K and Ka KQ75 masked maps, which is consis- 
tent with the fact that even these masked maps are still 
foreground contaminated at some level. 

The structure of the paper is as follows. In Sec. [Ill 
we introduce our non-Gaussianity indicators. Section Hill 
contains the results of applying our statistical indicators 
to the three-year and five-year WMAP data, and finally 
in Sec. [TVl we present the summary of our main results 
and conclusions. 



II. NON-GAUSSIANITY INDICATORS AND 
MAPS 

In this section we construct two statistical indicators 
that measure the large-angle deviation from Gaussianity 
of CMB temperature fiuctuation patterns. 

The main underlying idea in the construction of our 
non-Gaussianity indicators and the associated maps, is 
that the simplest ways of describing the deviation from 
symmetry about the CMB mean temperature, and a 
non-Gaussian degree of peakness are by calculating, re- 
spectively, the skewness S = /is/cr^, and the kurtosis 
K = /i4/(j^— 3 from the data, where /is and /i4 are 
the third and fourth central moments of the temperature 
anisotropics distribution, and a is the variance. Clearly 
calculating S and K for the whole celestial sphere would 
simply yield two dimensionless numbers describing the 
asymmetry of the data about the CMB average temper- 
ature. 

However, one can go a step further and obtain direc- 
tional information about deviation from Gaussianity if 
instead we take a discrete set of points {j = 1, . . . , A^d 
homogeneously distributed on the celestial sphere S'^ as 
the center of spherical caps of a given aperture 7 and cal- 
culate Sj and Kj for each cap. The values Sj and Kj can 
then be viewed as measures of the non-Gaussianity in the 
direction (6>j, of the center of the cap j . Such study 
of the individual caps can thus provide information {2Nc 
numbers) about possible large-angle violation of Gaus- 
sianity in the CMB data. A more systematic study can 
be made by taking the above set of points {j} as the 
center of the pixels for a homogeneous pixelization of 5*^, 
and by choosing caps of large-angle aperture to scan the 
whole celestial sphere with steps equal to the separation 



between the centers of adjacent pixels.^ In this way, we 
construct two scalar discrete functions S and K defined 
over the whole celestial sphere that encode measures of 
non-Gaussianity in CMB data. 

This constructive process can be formalized as fol- 
lows. Let Qj = Q{Oj^ <pj]j) G 5^ be a spherical cap, 
with an aperture of 7 degrees, centered at {Oj^(j)j)^ for 
j = 1, . . . , A^c- Define the scalar functions S : Qj M 
and K : Qj R, that assign to the j*^ cap, centered at 
{Oj^ 0j), two real numbers Sj and Kj given by 

^ J 2 = 1 

- ]vr^E(r^-r)'-3, (2) 

where is the number of pixels in the j*^ cap, Ti is 
the temperature at the i^^ pixel, T is the CMB mean 
temperature, and (jj is the standard deviation for each j. 

We now use the above homogeneously distributed 
points on S'^ to scan the celestial sphere with evenly dis- 
tributed spherical caps (of a chosen aperture 7) to cal- 
culate Sj and Kj for each cap. Clearly, the numbers Sj 
and Kj obtained in this way for each cap can then be 
viewed as a measure of non-Gaussianity in the direction 
of the center of that cap {Oj , ) . Patching together the 
Sj and Kj values for each cap, we obtain the indicators 
that are discrete functions S = S{0^ (j)) and K = K{0, (j)) 
defined over the celestial sphere, which measure the de- 
viation from Gaussianity as a function of direction (^, (j)). 
In this way, S = S{6,(j)) and K = K{0^(j)) give a scalar 
directional measure of non-Gaussianity over the celestial 
sphere. 

Now, since S = S{6, (j)) and K = K{6, (j)) are discrete 
scalar functions defined on S'^ they can also be viewed 
as maps of non-Gaussianity, and we can expand each of 
these functions in their spherical harmonics and calculate 
their angular power spectrum. Thus, for the skewness 
function S = S{0^(j))^ for example, one has 

00 I 

S{e,C^)=Y. bimY,m{0,cP) , (3) 

and can calculate the corresponding angular power spec- 
trum 

^^ = ^El^^™l'' (4) 

m 

in order to further quantify the angular scale informa- 
tion regarding the deviation from Gaussianity of CMB 



Note that here this pixehzation is only a practical way of choosing 
the centers of the caps homogeneously distributed on S"^. It is 
not related to the pixelization of the CMB maps. 
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SKEWNESS map ILC - Syr 




FIG. 2: Kurtosis indicator maps from the WMAP three (left 



data. Clearly, one can similarly expand the kurtosis 
function K = K{0^(j)) and calculate its angular power 
spectrum K^. It then follows that, if a large-scale non- 
Gaussianity is present in the original temperature distri- 
bution, it should significantly affect the S and K maps 
on the corresponding angular scales. 

In the next section we shall apply the indicators S = 
S{e, (j)) and K = K{e, (j)) to both WMAP three and five- 
year data. 



III. NON-GAUSSIANITY INDICATORS AND 
WMAP DATA 

The WMAP team have produced high angular resolu- 
tion maps of the CMB temperature fluctuations in five 
frequency bands: K-band (22.8 GHz), Ka-band (33.0 
GHz), Q-band (40.7 GHz), V-band (60.8 GHz), and 
W-band (93.5 GHz). Following the WMAP team rec- 
ommendation for tests of Gaussianity [l3], we used 
the Internal Linear Combination (ILC) maps of both the 
three-year and the five-year CMB data [l5| along with 
the new mask KQ75^ which is slightly more conservative 
than the KpO, i.e., KQ75 sky cuts is 28.4% while KpO 



SKEWNESS map ILC - Syr 




panel) and five-year (right panel) ILC maps with mask KQ75. 

KURTOSIS map ILC - Syr 




panel) and five-year (right panel) ILC maps with mask KQ75. 



cuts is 24.5% . We also test for non- Gaussianity the ILC 
full-sky five-year map along with the five frequency fore- 
ground uncleaned five-year maps with and without the 
KQ75 mask. In all cases we chose the HEALPix param- 
eter ^gi(^g = 256 which corresponds to a partition 
of the celestial sphere into 786 432 pixels. 

In our calculations of skewness and kurtosis indica- 
tor maps (hereafter referred to as ^S— map and i^— map) 
from each CMB map we have scanned the celestial sphere 
with spherical caps of aperture 7 = 90°, centered at 
A^c = 768, 3 072 and 12 288 points homogeneously gen- 
erated on the sphere by using HEALPix. However, to 
avoid repetition we only present a detailed analysis for 
Nc = 12 288 in the following. 

Figures [Hand [2] show, respectively, the Mollweide pro- 
jection of the S*— map and i^T— map in galactic coordi- 
nates obtained from the ILC WMAP three (left panels) 
and five-year (right panels) maps with the KQ75 mask. 
They clearly show that the S{0^ (j)) and K{0^ (j)) distri- 
butions of hot and cold spots (higher and lower values) 
for the indicators are not evenly distributed in the celes- 
tial sphere, suggesting at first sight non- Gaussianity of 
the ILC masked data. The comparison between the two 
6'— maps (Fig. [1]) and the two iC— maps (Fig. [2|) shows a 
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FIG. 3: Depicted the dipole, quadrupole and the remaining low I components for the skewness (left panels) and kurtosis (right 
panels) maps obtained from the five-year WMAP ILC with mask KQ75. 



great number of similarities for each pair of maps of the 
indicators, which is a very first indication of the robust- 
ness of our results with respect to the three and five-year 
WMAP data.^ Figures [T] and [2] are also suggestive of 
large-scale components in the maps of both indicators. 

To provide some additional qualitative information 
about the anisotropic distribution of our non-Gaussianity 



^ It is interesting to note the presence of great circles of unknown 
origin near the galactic plane in the K maps, and in a less ac- 
centuated way in the S maps. 



indicators, we depict in Fig. [3] the dipole and the 
quadrupole as well as the full S and K maps with these 
two components and the monopole removed. The left 
panels display these components for the skewness indica- 
tor S{0^(j)), while the right panels show the same com- 
ponents for the kurtosis indicator K{0^(j)). These mul- 
tipoles maps were calculated from the WMAP five-year 
ILG masked data, but corresponding maps for three-year 
WMAP data are largely similar to the depicted maps, 
and were not included to avoid repetition. 

It is known that the frequency K-band, Ka-band, Q- 
band, V-band, and W-band foreground uncleaned maps. 
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FIG. 4: Histograms show the distribution of the values for the multipoles 5*5 and calculated from the MC (scrambled) ILC 
maps. This figure illustrates the non-Gaussianity of a typical distribution for a given fixed i. 



have different contaminants. These features may appear 
in these maps in the form of non-Gaussianity even when 
the KQ75 mask is utilized. Thus, in order to suitably test 
the S and K indicators for non-Gaussianity one should 
calculate their maps and angular power spectra not only 
for the foreground reduced ILC masked and unmasked 
maps, but also for these single frequency maps with and 
without KQ75 mask. This also allows a comparative 
analysis of the outcomes. To this end, we have also calcu- 
lated the iS— maps and iC— maps for each of these five fre- 
quency maps with and without KQ75 mask, for caps with 
aperture 7 = 90°, centered at Nc = 12 288 points homo- 
geneously generated on the sphere by using HEALPix. 
However, to avoid repetition of figures which give only 
qualitative information, in the following we shall con- 
centrate on their angular power spectra, which provide 
quantitative information. 

We calculated the angular power spectrum of the S 
and K maps generated from the three and five-year data 
of ILC, and the foreground unreduced K, Ka, Q, V and 
W maps with and without KQ75 mask. These power 
spectra allow to estimate the statistical significance of 
Si and Ki by comparing them with the mean angular 
power spectrum of the S and K maps obtained from 
1 000 Monte- Carlo-generated (MC) statistically Gaussian 
CMB maps.'^ To make easier this comparison, instead of 
using the angular power spectra Si and K£ themselves, 
we employed the differential power spectra \S£ — Sil and 
\K£ — Ki\. Throughout the paper the mean quantities 
are denoted by over line. 

To describe with some details our calculations we fo- 
cus on the skewness indicator since similar proce- 



^ We note that each MC map is a stochastic realization of the 
WMAP best-fitting angular power spectrum of the ACDM 
model, obtained by randomizing the multipole temperature com- 
ponents airn within the cosmic variance limits 0, fisj . 



dure holds for the kurtosis indicator K. Starting from 
a given CMB seed map (ILC or any frequency band 
map) we generated 1 000 MC Gaussian {scrambled) CMB 
maps, which are then used to generate 1 000 skewness 
6'— maps, from which we calculate 1 000 power spectra: 
S\ with enumeration index i = 1, ••• ,1000. In this 
way, for each fixed multipole component Sl^^^^^ we have 
1 000 values from which we calculate the mean value 
'Si = (1/1000) 5']. From this MC process we 

have at the end ten mean values Si each of which are 
then used to compare with the corresponding power spec- 
trum component S^ of the skewness map obtained from 
an unscrambled seed temperature WMAP map, in order 
to evaluate the statistical significance of each multipole 
component. Thus, for example, to study the statisti- 
cal significance of the dipole moment of the skewness- 
ILC-map Sl^^ we calculate \Sl^^ — Si\^ where the mean 
dipole value Si is calculated from the i = 1, • • • , 1 000 
power spectra of the Gaussian {scrambled) maps. 

The panels of Fig. 2] display the histograms which show 
the distribution of the values for S^ and calculated 
from the scrambled ILC maps. This figure makes clear 
that a typical distribution of MC values for a given fixed 
£ is highly non-Gaussian. 

Figure [5] shows the differential power spectra calcu- 
lated from full-sky CMB five-year maps, i.e., it displays 
the absolute value of the deviations from the mean an- 
gular power spectrum of the skewness S^ (left panel) and 
kurtosis (right panel) indicators for ^ = 1, • • • , 10 , 
which is a range of multipoles values useful to investi- 
gate the large-scale angular characteristics of the S and 
K maps. This figure makes apparent the strong devia- 
tion from Gaussianity of the unmasked frequency maps. 
This is expected from the very outset since these full-sky 
band maps are highly contaminated. Figure [5] also shows 
a significant deviation from Gaussianity in five-year ILC 
unmasked data, i.e. the deviations \Si — S(^ \ and \K£ — Ki\ 
for the five-year ILC unmasked data are not within 95% 
of the MC value. Actually the values of Si and ob- 
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FIG. 5: Differential power spectrum of skewness \St — St \ and \Ki — Ki\ (left) and kurtosis (right) indicators calculated from the 
and five-year (right panels) WMAP CMB maps with no mask. The 68% and 95% confidence levels are indicated, respectively, 
by the dashed and dash-dotted lines. 



tained from the data are far beyond 95% off) the 
mean MC values (see Fig. [5]). These results quantify and 
make clear the suitability of the WMAP team recom- 
mendation to employ the new mask KQ75 for tests of 
Gaussianity. These spectra maps were calculated from 
the WMAP five-year data, but the corresponding spec- 
tra for three-year data are very similar to the depicted 
spectra. We stress that for the frequency maps we have 
used the foreground uncleaned maps as the temperature 
fluctuations seed maps. 

Figure [6] shows similar differential power spectra but 
now calculated from the CMB three and five-year maps 
with KQ75 mask. Tables H] and [II] complement Fig. [6] 
by collecting together the percentage of the deviations 
IS} — Si\ (calculated from 1 000 scrambled MC simulated 
maps) which are smaller than \S£ — Si] obtained from the 
data, i.e. from the five-year ILC, K, Ka, Q, V, and W 
masked maps. Thus, for example, according to Table [J 
for the K band 99.6% of multipoles obtained from 
the MC maps are closer to the mean S'3 than the value 
5*3 calculated from the K map (three and five-year data). 
This indicates how unlikely (only 0.4%) is the occurrence 
of the value obtained from the data for the multipole Ss 
in the set of MC simulated maps for this band, giving 
therefore a clear indication of deviation from Gaussianity 
for the K masked maps to the extent that this is not 
within - 95% of MC values. 

Regarding the angular power spectrum of the skew- 
ness and kurtosis maps calculated from the ILC maps 
with KQ75 mask, it is clear from the Fig. [6] and Tables [J 
and nil that the low £ multipoles (£=1, •••,10) values 

— indicative of large-angle deviation from Gaussianity 

— are not statistically significant, i.e. they are within 
the 95% of the MC values. In other words, the occur- 
rence of these ILC multipole values is very likely (> 5%) 



in the set of MC values.^ It is interesting to note that 
a comparison between Fig. [5] and Fig. [6] makes apparent 
the important role of the mask in the test for Gaussian- 
ity carried out. Furthermore, Fig. [6] and Tables [J and [III 
also show that this analysis is robust with respect to the 
three and five-year ILC maps. 

As regards the angular power spectra of S and K cal- 
culated from the K-band map with KQ75 mask, it is 
clear from Fig. [6] and Tables [J and [III that both the 
skewness and kurtosis spectra reveal deviations from 
the mean power spectrum values, greater than 95%, for 
£ = 2,3,5,6,7,8,10. This indication of deviation from 
Gaussianity comes chiefly from the fact that even with 
KQ75 masked K-band map is still foreground contam- 
inated as illustrated in first row [panels (a) and (b)] of 
Fig.lZl 

Figure [3 [second row, panels (c) and (d)] suggests a 
residual foreground contamination in the Ka-band map 
with KQ75 mask. This remnant contamination from 
galactic foreground emission is detected by the kurto- 
sis indicator whose power spectrum shows an excess 
of power relative to the mean for ^ = 1, 7, 9 components 
[see the second row of panels in Fig. [Hand Table HI] . 

Regarding the angular power spectra of the S and K 
indicators for the remaining three and five-year frequency 
maps (Q, V, W) with KQ75 mask the panels of Fig. [6] 
along with Tables [J and [11] show that the occurrence of 
the power spectrum low £ values calculated from these 
frequency maps are very likely (> 5%) in the set of MC 
values, i.e. the Q, V and W multipoles values are well 
within ^ 95% of MC values. This makes clear that for 
these frequency maps the KQ75 mask reduces consider- 
ably the contamination to a level which make the power 



^ We note that the values of Sq (five-year) and (five-year) which 
are slightly smaller than 5%. 
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FIG. 6: Differential power spectrum of skewness (first row) and kurtosis (second row) indicators calculated from the three (left 
panels) and five-year (right panels) WMAP CMB maps with mask KQ75. The 68% and 95% confidence levels are indicated, 
respectively, by the dashed and dot-dashed lines. See the text for more details. 
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K 




Ka 




Q 




V 




W 




ILC 


1 


90.6 — 


93.2 


83.7 — 


76.4 


88.5 — 


88.7 


74.5 — 


80.6 


76.0 — 


66.4 


14.8 — 


26.3 


2 


95.9 — 


97.7 


92.3 — 


90.1 


6.5 — 


1.3 


70.5 — 


22.8 


57.2 — 


23.2 


59.0 — 


57.9 


3 


99.6 — 


99.6 


24.1 — 


63.6 


10.0 — 


- 5.7 


66.1 — 


61.5 


57.7 — 


57.5 


96.3 — 


94.3 


4 


71.0 — 


94.4 


23.7 — 


16.0 


39.0 — 


34.3 


47.0 — 


46.6 


52.3 — 


45.6 


90.0 — 


83.7 
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96.4 — 


98.4 


6.4 — 


16.1 


38.8 — 


21.9 


74.6 — 


70.1 


80.8 — 


72.4 


86.0 — 


88.7 
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99.5 — 


99.7 


16.7 — 


30.1 


34.2 — 


23.6 


45.8 — 


42.5 


51.5 — 


43.3 


96.2 — 


96.0 
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99.7 — 


99.7 


24.4 — 


64.1 


53.0 — 


46.7 


88.7 — 


89.3 


87.5 — 


88.6 


33.4 — 


45.9 
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99.5 — 


99.7 


19.0 — 


16.0 


43.8 — 


43.4 


59.4 — 


57.7 


58.8 — 


57.8 


83.7 — 


59.4 


9 


99.6 — 


99.7 


41.6 — 


70.7 


53.4 — 


45.4 


87.1 — 


87.9 


87.2 — 


88.2 


47.0 — 


45.0 


10 


99.8 — 


99.8 


10.1 — 


- 4.1 


67.5 — 


63.4 


75.4 — 


77.6 


82.9 — 


81.9 


66.0 — 


72.3 



TABLE I: Percentage of the deviations \S\ — Si\ (for ^ = 1, • • • ,10 and calculated from 1 000 scrambled MC simulated maps), 
which are smaller than \Si — Si\ obtained from the data: 5*— maps generated from the three (left entries of each column) and 
five-year (right accesses for each column) CMB maps in the frequencies K (22.8 GHz), Ka (33.0 GHz), Q (40.7 GHz), V (60.8 
GHz), and W (93.5 GHz), along with the ILC three and five-year maps. 



spectra of this masked map compatible with Gaussianity. 

These results of our statistical analyses indicate that 
the current CMB temperature fluctuations of the ILC 
foreground reduced three and five-year KQ75 masked 
maps are consistent with Gaussianity in that the oc- 
currence of these ILC multipole multipole values is very 



likely (> 5%) in the set of MC values, i.e. the ILC mul- 
tipole values of the 5* and K are within :^ 95% of the 
MC values of the scrambled (Gaussian) maps. This result 
agrees with the WMAP team and other analyses made by 
using different statistical tools 0, B Hi • On the other 
hand, our analyses detect a deviation from Gaussianity 
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I 


K 




Ka 




Q 




V 




W 




ILC 


1 


99.9 — 


99.9 


98.7 — 


99.3 


44.5 — 


11.7 


82.0 — 


87.8 


86.4 — 


82.6 


19.4 — 18.9 


2 


99.9 — 


99.9 


27.1 — 


14.1 


29.8 — 


25.2 


38.4 — 


42.1 


41.6 — 


30.3 


74.9 — 76.4 


3 


99.9 — 


99.9 


1.7 — 


51.5 


90.1 — 


89.5 


86.9 — 


90.9 


88.7 — 


90.8 


16.0 — 10.9 


4 


99.8 — 


99.8 


3.0 — 


24.8 


55.4 — 


45.6 


84.8 — 


72.3 


89.1 — 


79.6 


37.2 — 42.9 


5 


99.9 — 


99.9 


7.9 — 


10.5 


69.4 — 


54.3 


80.1 — 


74.1 


93.2 — 


87.6 


94.7 — 95.3 


6 


99.8 — 


99.8 


21.0 - 


- 2.2 


66.1 — 


49.2 


73.0 — 


71.7 


80.1 — 


67.4 


8.3 — 1.9 


7 


99.9 — 


99.9 


95.1 — 


97.9 


86.5 — 


74.9 


92.5 — 


91.8 


93.7 — 


93.7 


60.3 — 53.1 


8 


99.8 — 


99.8 


7.2 — 


4.8 


82.8 — 


76.9 


91.8 — 


90.0 


93.2 — 


90.7 


4.3 — 13.4 


9 


99.9 — 


99.9 


94.8 — 


97.9 


5.9 — 


47.4 


85.6 — 


68.8 


92.3 — 


81.0 


93.5 — 95.0 


10 


99.8 — 


99.8 


3.2 — 


17.1 


48.7 — 


37.1 


89.2 — 


73.3 


80.4 — 


67.1 


2.2 — 2.3 



TABLE II: Percentage of the deviations \K\ — Ki\ (calculated from 1000 scrambled MC random maps for ^ = 1, • • • , 10), 
which are smaller than \Ki — Ki\ obtained from the data: K—maps generated from the three (left entries of each column) and 
five-year (right entries for each column) CMB seed maps in the channels K (22.8 GHz), Ka (33.0 GHz), Q (40.7 GHz), V (60.8 
GHz), and W (93.5 GHz), along with the ILC three and five-year maps. 




FIG. 7: CMB temperature fluctuations five-year maps in the frequencies K (22.8 GHz) and Ka (33.0 GHz) without (left panel) 
and with the mask KQ75. These figure illustrated that even these masked maps are still foreground contaminated at some 
level. These signs of deviation from Gaussianity in the KQ75 masked maps are captured by our indicators S (K KQ75 masked 
map) and K (K and Ka KQ75 masked maps). 



in the three and five-year frequency KQ75 masked K and 
Ka maps which is consistent with the fact that even these 
masked maps are foreground contaminated at some level 
(see Fig. [7]). As for the Q, V, and W frequency unre- 
moved foreground maps, our analyses indicate that al- 
though considerably contaminated in their full-sky form 
(as captured in Fig. [5]), KQ75 mask cuts the galactic re- 
gions so as to bring their S and K power spectra to a 



level consistent with Gaussianity (the multipoles values 
are within ~ 95% of MC values). 

Finally, to have an overall assessment of the power 
spectra of the S and K maps calculated from each CMB 
seed map we have performed a ^^st to find out the 
goodness of fit for Si and Ki multipole values as com- 
pared to the expected multipoles values from the Gaus- 
sian MC maps. In this way we can obtain one number 
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that collectively quantifies the extent to which a given 
KQ75 masked map is consistent with Gaussianity. For 
the power spectra we found that the ratio / dof 
(dof stands for degree of freedom) from the K, Ka, Q, 
V, W and ILC maps are given, respectively, by 21.5, 4.9, 
6.0, 5.2, 3.9 and 1.2, while for the kurtosis power spectra 
Ki of these maps the values of x^/dof are, respectively, 
35 652, 135, 0.5, 6.4, 5.6 and 0.4. Clearly the greater are 
these values the smaller are the probability, i.e., the 
probability that the multipoles values Si and (from 
each CMB maps) and the expected MC multipole val- 
ues agree. A cut-off value of x^/dof 7 to reject the 
the hypothesis that the multipole values Si and Ki are a 
good approximation to the expected multipole values of 
MC Gaussian realizations, is therefore enough to ensure 
that Q, V, W and ILC KQ75 masked are consistent with 
Gaussianity, while the K and Ka {KQ75 masked) maps 
are not. 

The calculations of our non- Gaussianity indicators re- 
quire not only the choice of a CMB map as input, but also 
the specification of some quantities whose choice could in 
principle affect the outcome of our results. To test the 
robustness of our scheme, hence of our results, we stud- 
ied the effects of changing various parameters employed 
in the calculation of our indicators. We found that the S 
and K angular power spectra do not change appreciably 
as we change the WMAP ILC three- year and five-year 
maps; the resolution of CMB temperature maps used 
(786 432 or 3145 728 pixels); and the number of point- 
centers of the caps with values 768, 3 072 and 12 288. 
The deviations from the mean power spectrum for three- 
year and five-year frequency bands maps are also robust 
(see Table |I] and Table ^ . 

Concerning the robustness of the above analyses some 
additional words of clarification are in order here. First, 
we note that the calculations of the S*— maps and 
iC— maps by scanning the CMB maps sometimes include 
caps whose center is within or close to the KQ75 masked 
region. In these cases, the calculations of these indicators 
are made with smaller number of pixels, which clearly 
introduce additional statistical noise as compared to the 
cases whose caps centers are far away from the mask. In 
order to minimize this effect we have scanned the CMB 
sky with spherical caps of aperture 7 = 90°.^ 



^ Note that even in the case of 90° caps the sky regions whose 
center are within or close to the galaxy have smaller sky fraction 
than those centered far away from the mask. This could lead 
to bias in our non-Gaussianity estimators. A possible way to 
circumvent this problem is by dividing Sj and Kj by the corre- 
sponding standard deviation {Sj — Sj'^Y^'^ or (K| - K^^)^/^, 
obtained by using MC CMB seed maps. We have calculated 
the maps and power spectra for these normalized indicators for 
1 000 MC maps. It turns out that for 7 = 90° caps, the power 
spectra of normalized and unnormalized indicators, calculated 
for the WMAP ILC five-year map, essentially coincide. For caps 
of smaller aperture as, e.g., 7 = 60°, although the normalization 
gives rise to smaller powers for and K^, the loss of powers for 



Second, we have checked the sensitivity of our results 
regarding ILC maps with respect to the angular resolu- 
tion of CMB maps, the number of point-centers of 
the caps used to scan the CMB sky, and the aperture 7 
of the caps, by carrying out correlation analyses between 
the resultant S and K maps calculated for CMB maps 
with 3145 728 and 786 432 pixels along with = 768, 
3 072 and 12 288. We found that for 7 = 90° the re- 
sulting S and K maps are strongly correlated according 
to Pearson's correlation coefficient. Thus, for example, 
for the S*— maps calculated from the two CMB pixeliza- 
tions the Pearson's correlation coefficients are equal to 
99.995%, 99.995% and 99.996% for A^c = 768, 3 072 and 
12 288, respectively, while for the i^— maps this coeffi- 
cient is 99.996% for the three values of Nc. However, for 
caps with aperture 7 = 70°, for example, the Pearson's 
correlation coefficients for S and K maps (calculated for 
CMB maps with these two pixelizations and the same 
values of Nc) fall below 50%, making clear that caps of 
aperture 7 = 90° are the most suitable spherical caps to 
deal with the statistical noise in our scheme. 



IV. FINAL REMARKS AND CONCLUSIONS 

We have proposed two new indicators for measuring 
large-angle directional deviations from Gaussianity in the 
CMB data, and have used them to search for the large- 
angle (low £) deviation from Gaussianity in the three and 
five-year foreground reduced ILC and the five single fre- 
quency (K, Ka, Q, V, W) KQ75 masked maps. Our 
directional indicators enable us to construct skewness S 
and kurtosis K maps (as, e.g.. Fig. [1] and Fig. [2]), mak- 
ing it possible to examine the presence and significance of 
possible large-angle non-Gaussianity in the WMAP CMB 
temperature fluctuations maps with and without KQ75 
mask. 

To obtain a more quantitative measure of non- 
Gaussianity we have studied the low £ angular power 
spectrum of the S and K maps generated from the (three 
and five-year) WMAP ILC and frequency (unremoved 
foreground) maps with and without KQ75 mask. For 
the full-sky ILC maps we found deviation from Gaus- 
sianity for the low ^, while for the ILC masked maps we 
found that the low £(^ = 1, •••,10) components are not 
significantly different from corresponding components of 
the expected power spectrum calculated from S and K 
maps obtained from 1 000 Monte Carlo CMB maps gen- 
erated by considering the Gaussian random hypothesis 
based on the concordance model [8]. Actually, we have 
found that the values of the multipoles Si and Ki (for 
^ = 1, • • • , 10 ) are not statistically significant, i.e. they 
are within the 95% values of 5'^'s and Ki's of the MC 



each i is not enough to make the 7 = 60° results close to the 
clearly statistically less noisy 7 = 90° power spectrum estimates. 
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randomly scrambled maps] (Fig. [6] and Tables HI and . 
As regards the frequency maps, we found clear indica- 
tion of deviation from Gaussianity in the three and five- 
year frequency KQ75 masked maps: K and Ka, which 
is expected and consistent with the fact that even these 
masked maps present some level of foreground contam- 
ination away from the masked region (see Fig. [7|). The 
deviation for the Q, V, and W masked maps are within 
the 95% expected values from MC randomly simulated 
maps. 

To have an overall quantitative assessment of the power 
spectra Si and calculated from K, Ka, Q, V, W and 
ILC maps we have performed a t^st to determine the 
goodness of fit for low £ multipole values as compared 
to the expected multipoles values from the Gaussian MC 
maps. In this way we have obtained numbers that collec- 
tively quantify the extent to which these KQ75 masked 
maps are consistent with Gaussianity. 

The results of our statistical analyses indicate that the 
current CMB temperature fiuctuations ILC three and 
five-year masked data are consistent with Gaussianity, 
in agreement with the WMAP team and other analyses 
made by using different statistical tools [gl, l7|, |8|]. We 
have demonstrated that the results of our analyses are 
robust by showing that the iS— map and iC— map do not 
significantly change with different choices of variables in- 
volved in our scheme, so long as the statistical noise is 
kept under control. 

The effects of different foreground-reduced algorithms 
as detected by our non- Gaussianity indicators for other 
WMAP three and five-year maps [11, 0, |^ is un- 
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